Draft version February 1, 2008 

Preprint typeset using LAT^^ style cniulateapj v. 26/01/00 



RAPIDLY ROTATING STRANGE STARS FOR A NEW EQUATION OF STATE OF STRANGE 

QUARK MATTER 



Ignazio Bombaci"'^, Arun V. Thampan^ 



Bhaskar Datta^'" 



o 
o 
o 

o 



> 

oo 
(N 

m 

a\ 
o 
o 
o 

Oh- 
I 

o 



% 



Draft version February 1, 2008 

ABSTRACT 

For a new equation of state of strange quark matter, we construct equilibrium sequences of rapidly 
rotating strange stars in general relativity. The sequences are the normal and supramassive evolutionary 
sequences of constant rest mass. We also calculate equilibrium sequences for a constant value of f2 
corresponding to the most rapidly rotating pulsar PSR 1937 +21. In addition to this, we calculate 
the radius of the marginally stable orbit and its dependence on 17, relevant for modeling of kilo-Hertz 
quasi-periodic oscillations in X-ray binaries. 

Subject headings: stars: neutron - equation of state - general relativity 



1. INTRODUCTION 

One of the most fascinating aspects of modern astro- 
physics, is the possible existence of a new family of col- 
lapsed stars consisting completely of a deconfined mixture 
of up (m), down (d), and strange (s) quarks, together with 
an appropriate number of electrons to guarantee electrical 
neutrality. Such compact stars have been referred to in 
the literature, as strange stars (SS), and their constituent 
matter as strange quark matter (SQM). The possible ex- 
istence of SS is a direct consequence of the the so called 
strange matter hypothesis^ formulated by Witten (1984) 
(see also Bodmer 1971, Terazawa 1979). According to this 
hypothesis, SQM, in equilibrium with respect to the weak 
interactions, could be the absolute ground state of strongly 
interacting matter rather than ^^Fe. Ever since the formu- 
lation of this hypothesis, there have been several reports 
in the literature, on the structure of non-rotating SS (e.g. 
Haensel, Zdunik & Schaefer 1986; Alcock, Farhi & Olinto 
1986). However, these have remained largely theoretical 
and speculative in nature. Observational data from RXTE 
of certain low mass X-ray binaries (LMXBs) have recently 
given rise to a resurgence in the astrophysical interest as- 
sociated with SS. It has been suggested that at least a few 
LMXBs could be harboring SS. For example, recent stud- 
ies have shown that the compact objects associated with 
the X-ray burster 4U 1820 -30 (Bombaci 1997), the burst- 
ing X-ray pulsar GRO J1744 -28 (Cheng et al. 1998) and 
the X-ray pulsar Her X-1 (Dey et al. 1998) are likely SS 
candidates. The most promising and convincing SS candi- 
dates are the compact objects in the bursting millisecond 
X-ray pulsar SAX J1808.4 -3658 (Li et al. 1999a), and in 
the atoll source 4U 1728 -34 (Li et al. 1999b). 

The compact nature of these sources make general rel- 
ativity important in describing these systems. Further- 
more, their existence in binary systems imply that these 
may possess rapid rotation rates (Bhattacharya & van den 
Heuvel 1991 and references therein). Particularly, the two 
SS candidates in SAX J1808.4 -3658 and 4U 1728 -34 
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are millisecond pulsars having spin periods P = 2.49 ms 
and P = 2.75 ms respectively. These two properties make 
the incorporation of general relativistic effects of rotation 
imperative for satisfactory treatment of the problem. 

Most of the calculations on the rotational properties of 
SS, reported so far, have relied on the slow rotation ap- 
proximation (Colpi & Miller 1992, Glendenning & We- 
ber 1992). This approximation loses its validity as the 
star's spin frequency approaches the mass shedding limit. 
Rapidly rotating SS sequences have been recently com- 
puted by Gourgoulhon et al. (1999), Stergioulas et al. 
(1999), and Zdunik et al. (2000). However, all the calcula- 
tions mentioned above make use of a very schematic model 
(Freedman & McLerran 1978, Farhi & Jaffe, 1984), related 
to the MIT bag model (Chodos et al. 1974) for hadrons, 
for the equation of state (EOS) of SQM. Within the MIT 
bag model EOS, the SS radii calculated are seen to be in- 
compatible with the mass-radius {M~R) relation (Li et al. 
1999a) for SAX J1808.4 -3658, and only marginally com- 
patible (see Fig. 2) with that for 4U 1728 -34 (Li et al. 
1999b). 

In this letter, we present calculations of equilibrium se- 
quences of rapidly rotating SS in general relativity using 
a new model for the EOS of SQM derived by Dey et al. 
(1998). This model is based on a "dynamical" density- 
dependent approach to confinement. In contrast, in the 
simple EOS based on the MIT bag model, the medium 
effects on the quark degrees of freedom and on the quark- 
quark interaction are not considered. For illustrative pur- 
poses, we compare our results, obtained using this new 
EOS, with those for the MIT bag model EOS. 

We use the methodology described in detail in Datta, 
Thampan & Bombaci (1998) to calculate the structure of 
rapidly rotating SS. For completeness, we briefly describe 
the method here. For a general axisymmetric and station- 
ary space-time, assuming a perfect fluid configuration, the 
Einstein field equations reduce to ordinary integrals (using 
Green's function approach). These integrals may be self 
consistently (numerically and iteratively) solved to yield 

via Buonarroti 2, 1-56127, Italy; BOMBACI@pisa.infn.it 
Pune 411 007, India; arunOiucaa.ernet.in 



the value of metric coefficients in all space. Using these 
metric coefficients, one may then compute the structure 
parameters, angular momentum and moment of inertia 
corresponding to initially assumed central density and po- 
lar to equatorial radius ratio. These may then be used 
(as described in Thampan & Datta 1998) to calculate pa- 
rameters connected with stable circular orbits (like the in- 
nermost stable orbit and the Keplerian angular velocities) 
around the configuration in question. 

The sequences that we calculate are: constant rest mass 
sequences, constant angular velocity sequences, constant 
central density sequences and constant angular momen- 
tum sequences. We also calculate the radius rorb of the 
marginally stable orbit and its dependence on the spin 
rate of the SS, which will be relevant for modeling X-ray 
burst sources involving SS. 

2. THE EQUATION OF STATE FOR SQM 

As mentioned earlier, the schematic EOS for SQM based 
on the MIT bag model, has become the "standard" EOS 
model for SS studies. However, this EOS model becomes 
progressively less trustworthy as one goes from very high 
density region (asymptotic freedom regime) to that of 
low densities, where confinement (hadron formation) takes 
place. Recently, Dey et al. (1998) derived a new EOS for 
SQM using a "dynamical" density-dependent approach to 
confinement. This EOS has asymptotic freedom built in, 
shows confinement at zero baryon density, deconfinement 
at high density. In this model, the quark interaction is de- 
scribed by a colour-Debye-screened inter-quark vector po- 
tential originating from gluon exchange, and by a density- 
dependent scalar potential which restores chiral symmetry 
at high density (in the limit of massless quarks). This 
density-dependent scalar potential arises from the density 
dependence of the in-medium effective quark masses Mq, 
which are taken to depend on the baryon number density 
n-Q according to (see Dey et al. 1998) 
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where no — 0.16 fm^'^ is the normal nuclear matter den- 
sity, q(= u, d, s) is the flavor index, and j/ is a parameter. 
The effective quark mass Mq(nB) goes from its constituent 
mass at riB = 0, to its current mass rriq, as ne ^ cxo. Here, 
we consider two different parameterizations of this EOS, 
which correspond to a different choice for the parameter v. 
The equation of state SSI (SS2) corresponds to j^ = 0.333 
[y = 0.286). These two models for the EOS give absolutely 
stable SQM according to the strange matter hypothesis. 

In order to compare our results with those of previous 
studies, we also use the MIT bag model EOS, for massless 
non-interacting quarks and B = 90 McV/fm"^ (hereafter 
the B90_0 EOS). 

3. RESULTS AND DISCUSSION 

The equilibrium sequences of rotating SS depend on two 
parameters: the central density (pc) and the rotation rate 
($7) . For purpose of illustration, we choose three limits in 
this parameter space. These are: (i) the static or non- 
rotating limit, (ii) the limit at which instability to quasi- 
radial mode sets in and (iii) the centrifugal mass shed limit. 
The last limit corresponds to the maximum f2 for which 



centrifugal forces are able to balance the inward gravita- 
tional force. 

The result of our calculations for EOS SSI is displayed in 
Fig. 1. In Fig. 1 (a) we show the functional dependence of 
the gravitational mass (M) with pc- In these set of figures, 
the bold solid curve represents the non-rotating or static 
limit, and the bold dashed curve the centrifugal mass shed 
limit. The thin solid curves are the constant rest mass 
(Mo) evolutionary sequences. The evolutionary sequences 
above the maximum stable non-rotating mass configu- 
ration arc the supramassivc evolutionary sequences, and 
those that lie below this limit arc the normal evolution- 
ary sequences. The maximum mass sequence for this EOS 
corresponds to Mq = 2.2 Mq . The thin dot-dashed line 
(slanted towards left) represents instability to quasi-radial 
perturbations. In the central panel of the same figure, we 
give a plot of M as a function of the equatorial radius R. 

In panel (c) of Fig. 1, we display the plot of fJ as a 
function of the specific angular momentum j = cJ/GMq 
(where J is the angular momentum of the configuration). 
Unlike for neutron stars (e.g. Cook, Shapiro & Teukolsky 
1994; Datta, Thampan & Bombaci 1998), the fi-j curve 
does not show a turn-over to lower j values for SS. This is 
due to the effect of the long-range (non-perturbative) in- 
teraction in QCD (Quantum Chromo-Dynamics), which is 
responsible for quark confinement in hadrons, and makes 
low mass SS self-bound objects. 17 for the mass shed limit 
appears to asymptotically tend to a non-zero value for 
rapidly rotating low mass stars. A further ramification of 
this result is that the ratio of the rotational energy to the 
total gravitational energy (T/W) becomes greater than 
0.21 (as also reported by Gourgoulhon et al. 1999) thus 
probably making the configurations susceptible to triaxial 
instabilities. 

In Table 1 we display the values of the structure pa- 
rameters for the maximum mass non-rotating SS models. 
The larger value of the maximum mass for the SSI model, 
with respect to the SS2 model, can be traced back to role 
of the parameter ly in eq. (1) for the effective quark mass 
Mg. In fact, a larger value oi v (SSI model) gives a faster 
decrease of Mq with density, producing a stiffer EOS. 

Table 2 and Table 3 displays the maximum mass ro- 
tating and maximum angular momentum models for the 
EOS models under consideration. While for EOS SSI, the 
maximum mass rotating model and the maximum angular 
momentum models are the same, for EOS SS2, the two 
models are slightly different, with the maximum angular 
momentum model coming earlier (with respect to pc) than 
the maximum mass rotating configuration. 

In Table 4 we list the values of the various parameters 
for the constant fl sequences for EOS SSI. The first entry 
in this table corresponds to pc for which rorb = R- For 
higher values of pc, rorb > Rj for large values of pc, the 
boundary layer (separation between the surface of the SS 
and its innermost stable orbit) can be substantial (^ 5 km 
for the maximum value of the listed pc) . 

In Fig. 2, we plot our theoretically calculated M-R 
curves (dashed curves) for EOS SSI and B90_0, for a ro- 
tational frequency of 364 Hz corresponding to the inferred 
rotational frequency of the compact star in the source 
4U 1728 -34 (Mendez & Van der Khs 1999). In the same 
figure, we plot the radius Rq of the inner edge of the ac- 



cretion disk (dotted curve) for 4U 1728 —34 as deduced by 
Titarchuk & Osherovich (1999) (see also Li et al. 1999b) 
from a fit of kHz QPO data in this source. Since Rq must 
be larger than both R and Torb, it is possible to deduce (Li 
et al. 1999b) an upper bound for the radius and mass of 
the central accretor. Using these constraints on M and R, 
Li et al. (1999b) concluded that 4U 1728 -34 is possibly 
a strange star rather than a neutron star. However, in 
their calculation, Li et al. (1999b) used an approximate 
EOS-independent expression to account for the effects of 
rotation on the moment of inertia of the star and hence on 
?'orb (see Li et al 1999b for further details). In contrast, in 
this letter we make an "exact" calculation (as a result of 
our accurate calculation of angular momentum and hence 
the moment of inertia) of rorb (triple-dot dashed curves) 
for the two EOS SSI and B90_0. Notice that the upper 
bound for the mass of 4U 1728 —34 (the intersection point 
between the curves Rq and rorb) is about 1.15 Mq, in 
agreement with the results of Li et al. (1999b). The two 
EOS models considered in Fig. 2 are both consistent with 



4U 1728 —34. Also, having in mind the scaling with B ^/^ 
of the M-R relation of SS within the bag model EOS (Wit- 
ten, 1984) and the constraints on the allowed values of the 
constant B to fulfil the strange matter hypothesis (Farhi 
& Jaffe, 1984), we see that the simple EOS based on the 
MIT bag model is only marginally compatible with the 
M-R relation for 4U 1728 -34. 

To summarize, in this letter, we present calculations 
of equilibrium sequences of rapidly rotating strange stars 
in general relativity, for a set of new equations of state. 
We compare the results so obtained with those for the 
generally used MIT bag model EOS. In addition to this, 
we illustrate the results of our computations for a spe- 
cific rotation rate inferred for the central accretor in the 
source 4U 1728 —34. We find that if the compact star 
in this source were indeed to be a strange star, then its 
mass and radius would be bounded above by a value of 
M = 1.15 Mq and R = 9.3 km. It is expected that future 
observations may shed more light on this issue. 
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Table 1 
Structure parameters for the non-rotating maximum mass configurations. Listed, are the central density (pc 

IN units of G CM~^, the gravitational mass (M) IN SOLAR UNITS, THE EQUATORIAL RADIUS {R) IN KM, THE BARYONIC 
mass (Mo) IN SOLAR UNITS, THE RADIUS OF THE MARGINALLY STABLE ORBIT (rorb) IN KM AND THE MOMENT OF INERTIA (7) IN 

UNITS OF 10'' 



EOS 


Pc 


M 


R 


Mo 


^orb 


/ 


SSI 

SS2 

B90.0 


4.65E+15 
5.60E+15 
3.09E+15 


1.438 
1.324 
1.603 


7.093 
6.533 

8.745 


1.880 
1.658 
1.937 


12.729 
11.716 
14.020 


0.736 
0.578 
1.099 



Table 2 
Structure parameters for the maximally rotating (fl = nms) maximum mass configuration. In addition to the 

QUANTITIES LISTED IN THE PREVIOUS TABLE, WE DISPLAY, THE ROTATION RATE {Q.) IN 10'* RAD S"\ THE RATIO OF THE 
ROTATIONAL TO THE TOTAL GRAVITATIONAL ENERGY {T/W) AND THE SPECIFIC ANGULAR MOMENTUM (j = cJ/GM^). 

EOS Pc n I M T/W R j Torb Mq 



SSI 

SS2 

B90_0 


3.10E+15 
3.60E+15 
1.90E+15 


1.613 
1.738 
1.190 


2.072 
1.613 
3.369 


2.077 
1.904 
2.272 


0.219 
0.218 
0.232 


10.404 

9.612 

13.213 


0.524 
0.570 
0.633 


11.656 
10.758 
14.612 


2.694 
2.366 
2.683 



Table 3 
Structure parameters for the maximum angular momentum configuration. 



EOS 


Pc 


n 


/ 


M 


T/W 


R 


J 


rorh 


Mo 


SSI 

SS2 

B90_0 


3.10E+15 
3.40E+15 
1.70E+15 


1.613 
1.719 
1.161 


2.072 
1.633 
3.456 


2.077 
1.899 
2.254 


0.219 
0.220 
0.239 


10.404 
9.693 
13.447 


0.524 
0.575 
0.650 


11.656 
10.837 
14.864 


2.694 
2.355 
2.650 



Table 4 
Structure parameters for the constant angular velocity sequence for EOS SSI. This sequence corresponds to 

THE rotation RATE OF THE PULSAR PSR 1937 +21 (BACKER et al. 1982), HAVING f2 = 4.03 X 10^ RAD S"^ OR PERIOD 

P= 1.556 MS. 



M T/W R j Torb Mq 



1.70E+15 


0.353 


0.852 


0.013 


6.663 


0.153 


6.663 


1.027 


1.80E+15 


0.433 


0.963 


0.012 


6.884 


0.143 


7.664 


1.178 


1.90E+15 


0.502 


1.052 


0.019 


7.036 


0.136 


8.378 


1.301 


2.40E+15 


0.701 


1.297 


0.010 


7.326 


0.116 


10.355 


1.660 


2.60E+15 


0.737 


1.346 


0.010 


7.344 


0.112 


10.741 


1.735 


4.60E+15 


0.769 


1.458 


0.008 


7.139 


0.096 


11.732 


1.914 


5.65E+15 


0.734 


1.449 


0.007 


7.007 


0.093 


11.703 


1.899 



Fig. 1. — Structure parameters for rotating strange stars corresponding to EOS SSI. The bold-solid line represents the non-rotating limit, 
the bold-dotted line the mass-shed limit and the almost vertical thin dot— dashed line tilted to the left is the instability limit to quasi-radial 
mode perturbations. The thin solid lines (labelled 1, 2 ...) represent constant rest mass sequences: 1: 1.59 Mq , 2: 1.66 Mq , 3: 1.88 Mq , 
4: 2.14 Mq 5: 2.41 Mq . 
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Fig. 2. — Mass— radius relations and r^^^, for a rotation rate of f2 = 2.287 X 10"^ rad s^^ (i.e. 364 Hz) corresponding to the inferred 
rotation rate of the central accretor in 4U 1728 —34 (Medez and Van der Klis 1999). Also plotted are the bounds on the M-R from general 
considerations, and the fit for these values obtained by Titarchuk & Osherovich (1999). See text for further details. 
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